Breathing mode frequencies of a rotating Fermi gas in the BCS-BEC crossover region 
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We study the breathing mode frequencies of a rotating Fermi gas trapped in a harmonic plus 
radial quartic potential. We find that as the radial anharmonicity increases, the lowest order radial 
mode frequency increases while the next lowest order radial mode frequency decreases. Then at a 
critical anharmonicity, these two modes merge and beyond this merge the cloud is unstable against 
the oscillations. The critical anharmonicity depends on both rotational frequency and the chemical 
potential. As a result of the large chemical potential in the BCS regime, even with a weak anhar- 
monicity the lowest order mode frequency increases with decreasing the attractive interaction. For 
large enough anharmonicities in the weak coupling BCS limit, we find that the excitation of the 
breathing mode frequencies make the atomic cloud unstable. 
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I. INTRODUCTION 

The rapid progress of ultra-cold atomic gas experi- 
ments provides unique opportunities for well controlled 
studies of quantum many body physics. For the case 
of Fermi atomic systems, the possibility of controlling 
the s-wave scattering length (a) between two different 
spin components allows to control the interaction by us- 
ing a magnetically tuned Feshbach resonance ,1]. This 
unique capability allows one to investigate the cross-over 
between the weakly interacting BCS regime (the regime 
where a — > CP) and Bose-Einstein condensate of dimers 
(the regime where a — ► + ) 2]. These two regimes meet 
in strongly interacting limit where the scattering length 
is divergent and at this unitarity limit, the physics is ex- 
pected to be universal Q. 

The appearance of quantized vortices of a quantum 
fluid under rotation offers direct evidence of superfluid- 
ity For example, the observation of quantized circula- 
tion in a rotating superfiuid 4 He Q and the observation 
of vortex lattice in a rotating Fermi gas of 6 Li [B| are 
two classic demonstrations of phase coherence in a su- 
perfiuid. These are analog to the vortex lattice in type- 
II superconductors in the presence of a magnetic filed. 
These vortices melt as the magnetic field increases and 
then the superconductors turn into normal at sufficiently 
large magnetic fields. For the case of rapidly rotating 
Fermi gasses, the force due to the trapping frequency al- 
most balances the centrifugal force and superfiuid cloud 
spreads in the plane perpendicular to the rotation axis. 
At the limit of very large rotation, the theory predicts 
that the atomic system enters into the fractional quan- 
tum Hall regime [g, 0|- However, the fractional quantum 
Hall window is expected to be very small and inversely 
proportional to the number of atoms in the trap. A possi- 
ble way of stabilizing the fractional quantum Hall regime 
is to add a positive quartic trapping potential. 

In this paper, we study the collective breathing mode 
frequencies of a rotating Fermi gas in the presence of a 
quartic trapping potential by using a hydrodynamic ap- 
proach. A negative, but small quartic term is always 



present with the Gaussian optical potentials in current 
experimental setups while added positive quartic term 
ensures the stability of the fast rotating regime. Ther- 
modynamic properties of a Bose gas confined in har- 
monic plus quartic potential trap can be found in ref. [|| . 
As breathing mode frequencies are very sensitive to the 
equation of state, these dynamical quantities can be used 
as tests for various theories. The breathing mode fre- 
quencies have been measured for non-rotating Fermi sys- 
tems in the BCS-BEC region [j| [l(| • In most of the pa- 
rameter regions, experimental results agree well with the 
hydrodynamic approaches, variational approaches and 
sum rule approaches However, the measured finite 



temperature axial and radial breathing mode frequen- 
cies show a striking increase in the intermediate BCS 
regime [tl, [l(| , in contrast to zero temperature theoreti- 
cal calculations in a harmonic trap. Furthermore, exper- 
imentalists were unable to measure the breathing mode 
frequencies in the weak coupling BCS limit. The devia- 
tion of the experimental data from theoretical results and 
the lack of experimental data in the weak coupling BCS 
limit were believed to be due the large Landau damping 
when the superfiuid energy gap is much smaller than the 
collective oscillation energies. With inclusion of a posi- 
tive quartic term in the trapping potential, we find some- 
what similar deviation of the breathing mode frequencies 
in the intermediate regions of BCS regime. We find that 
the breathing mode frequencies deviate significantly from 
the modes frequencies calculated in a harmonic trap and 
the atomic cloud is unstable against the breathing mode 
oscillations at larger chemical potentials. As the chemical 
potential is larger in the weak coupling BCS limit, exci- 
tation of breathing mode frequencies make the atomic 
cloud unstable. It should be noted that the Gaussian op- 
tical trap potential provides a negative quartic term in 
the trapping potential in experimental setups. We inves- 
tigate the effect of negative quartic term and find that 
the breathing mode frequencies tend to decrease in the 
entire BCS-BEC crossover region. 

This paper is organized as follows. In section II, we 
present the derivation of breathing mode frequencies us- 



2 



ing a hydrodynamic approach. In section III, we present 
our results together with a discussion. Finally in section 
IV, we draw our conclusions. 



II. FORMALISM 

We consider a rotating Fermi atomic system trapped in 
a harmonic plus radial quartic potential in the BCS-BEC 
crossover region. The trapping potential is 



V ex (r,z) = hlioy 



(1) 



where M is the atom mass, u>i's are the harmonic trap- 
ping frequencies, and r 2 = x 2 + y 2 . The Fermi atomic 
system rotates about the z-axis at frequency f2. 

We restrict ourselves to the case of large number of vor- 
tices in the system where the wavelength of the oscillation 
frequencies is much larger than the inter-vortex distance. 
This condition always satisfies as the typical wavelength 
of the lowest mode oscillations is of the order of the sys- 
tem size. Further, we assume that the vortices are uni- 
formly distributed in the superfluid so that we do not 
have to consider microscopic details of the single vortices. 
These assumptions are valid at the limit of large rotations 
where the atomic cloud spreads in the plane perpendicu- 
lar to the rotation. Within this diffused vorticity approx- 
imation 12j, diffuse vorticity is given by V x v = 2f2, 
where superfluid velocity is given by v = (H/2M)\79. 
The local superfluid density n and the local phase 9 are 
related through the wave function ip = ^/ne %e . The uni- 
form vortex density is given by n v = 2MQ/H. 

Assuming local equilibrium, we start with the continu- 
ity and Euler equations of rotational hydrodynamics, 



dn 
~dt 



-V • [n(r)v] 



(2) 



and 



71/ 



-V[-Mv 2 + V ex (r,z) 



-2Mv x fi + Mv x V x v (3) 



This hydrodynamic description is valid as long as the 
collisional relaxation time r is much smaller than the in- 
verse of the oscillation frequencies; cut « 1. The equa- 
tion of state enters through the density dependent local 
chemical potential (i(n). We fix the local chemical po- 
tential by introducing the equation of state in the form 
of fi{n) a n 1 . As we will discuss in the next subsec- 
tion, the polytropic index 7 is calculated by the method 
proposed by Manini and Salasnich [l3| in the entire BCS- 
BEC crossover region. Linearizing the density n and the 
superfluid velocity v around their equilibrium values as 
n = no(r) + Sn, v = vq + <5v and /i(n) = fx(no) + S/i with 



5fx = (dfi/dn)\ n — no Sn, we obtain the linearized version 
of the hydrodynamic equations. 



dSn 



-V • [no(r)*v] 



and 



2M5\ x n 



(4) 



(5) 



Starting from these two linearized equations, collec- 
tive breathing mode frequencies have been calculated in 
ref.[3] and ref.[l5| for a harmonic trap. As the authors 
have used two different ansatz for the velocity fluctua- 
tion, they produce two different results for the breathing 
mode frequencies. In this paper we closely follow the ap- 
proach adopted in ref . [l5| . generalizing the theory to an 
inharmonic trap. The^ate used in ref.0 ensures the 
conservation of angular momentum properly. In order 
to solve the linearized equations for the breathing mode 
frequencies, we take the equilibrium density in the lo- 
cal density approximation as n (r) cx [^0 — (l/2)-^(^ 2 — 
n 2 )r 2 - (1/2)Muj 2 z z 2 - ( K/4)r 4 } 1 '~< and use following vari- 
ational ansatz for the density fluctuations and velocity 
fluctuations. 



8w = {8Q.X x r + Sfl 2 x r 2 r + V[aj_r 2 + a z z 2 + f3r 4 ]} 

x exp[— iujt](6) 



and 



Sn = 



1-7 



{a 



a±r 



a z z 



br 4 } expHwi] (7) 



The first two terms 5£li and Stl? in Eq. ([6|) are parallel to 
the axis of rotation and guarantee that angular momen- 
tum is conserved during the oscillations. Substituting 
these two ansatz into Eq. and Eq. (O, we derive four 
linear equations for the variational parameters ao, a±, 
a 2 , and b. These linear equations yield three non-zero 
solutions for the breathing mode frequencies w m = uj/uj r 
as roots of the following equation. 



A + Buj 2 



Cut 



= 



(8) 



with, A = -(1 - C 2 ) 2 <5 2 [64d 7 ( 7 + 1) + 48 7 2 + 56 7 + 
16]-(l-C 2 )C 2 <5 2 (32 7 2 + 104 7 + 48)-16C 4 5 2 ( 7 + 2), B = 
(l-C 2 ) 2 [8( 7 +l)(2 7 + l) + 16d 7 ( 7 + 2)] + (l-C 2 )[<5 2 (8 7 2 + 
26 7 + 12) + C 2 (40 7 + 24)] + 8C 2 <5 2 ( 7 + 2) + 16( 4 and C = 
(l-C 2 )(2-10 7 )-5 2 ( 7 + 2)-8. The constants, C = 0/w r , 
S = u z /uj r , andd = [l/{l-C 2 ) 2 }(Kh/M 2 Lu?){n /huj r ) are 
a set of dimensionless parameters. The three solutions of 
Eq. (JSj) are the lowest order axial breathing mode fre- 
quency u}\ and the lowest and next lowest order radial 
breathing mode frequencies u>2 and uj%. 
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The Effective polytropic index and the chemical 
potential in the BCS-BEC crossover region 

We use the proposal made by Manini and Salasnich [l3| 
to calculate the effective polytropic index 7 and the 
chemical potential \x in the BCS-BEC crossover region. 
In the weak coupling BCS limit (a — > - ) and the uni- 
tarity limit (a — ► 00), the polytropic index is 7 = 2/3. 
In the deep BEC limit (a — ► + ), the polytropic index is 
7 = 1. In the BCS-BEC crossover regime, the scattering 
length's dependence on 7 is given by [l3| 



7 



2/3-2ye'(y)/5 + y 2 e"(y)/15 

e(y) - ye'(y)/5 



(9) 



with the parameter y — l/(kfd) is the interaction pa- 
rameter with kf being the Fermi wave vector. The func- 
tion e(y) is related to the energy per atom given by 
E = (3/5)E f e(y), where Ej = H 2 k 2 f /2M is the Fermi 
atomic energy of a non-interacting Fermi system in the 
trap. Above e'(y) = de(y)/dy and the double prime in- 
dicates the second derivative of the function on its argu- 
ment. Using the data presented in reference (l6| . Manini 
and Salasnich [l3| used a data fitting scheme to derive 
an analytical form of the function e(y) in the entire BCS- 
BEC region, 



e(y) = ol\ — OLi arctan 



ft 



\y\ 



(10) 



Two different sets of parameters are proposed for ccj's and 
ft's m the BCS regime (y < 0) and the BEC regime (y > 
0). In the BCS regime, the parameters are a\ = 0.42, 
a 2 = 0.3692, a 3 = 1.044, ft = 1.4328, and ft = 0.5523. 
In the BEC regime, the parameters are a.\ = 0.42, 012 = 
0.2674, a 3 = 5.04, ft = 0.1126, and ft = 0.4552. The 
expression for the chemical potential fi is given by [l3[ 



H = E f [e(y)-ye'(y)/5} 



(11) 



We determine the Fermi energy Ef of a non- interacting 
Fermi system in a harmonic plus radial quartic potential 
through the number equation. 




FIG. 1: Central chemical potential (fio) of a non-rotating 
Fermi g £LS ctS ct function of interaction parameter l/(fe/a) for 
various values of K. From top to bottom K — 0.05 (long black 
dashed line), and 0.01 (short black dashed line), harmonic 
trap (black solid line), —0.005 (short gray dashed line), and 
—0.01 (long gray dashed line). For the calculation, we use 
N = 2.0 x 10 6 and S = 0.045. 



2.05 
2.00 
1.95 
1.90 
1.85 
1.80 





/ 
J 









-1.5 



-1.0 



-0.5 



0.0 

l/(*/a) 



0.5 



1.0 



1.5 



FIG. 2: The lowest order radial breathing mode frequencies 
of a rotating Fermi gas in a harmonic trap with aspect ratio 
S = 0.045. The rotation frequencies are £ = (solid line), 0.5 
(long dashed line) and 0.9 (short dashed line). 



where the scaled parameters are K = hK/(M 2 u!^) and 
Ef = Ef/(huj r ). The upper and lower signs are cor- 
responding to K > and K < respectively. For 
the case of harmonic trap, the Fermi energy is Ef = 
^ r [37V^(l-C 2 )] 1/3 . 



N 



1 M 2 ojI ( E 



5/2 



f(Ef) 



15J y h\K\ \hu>r 

We define the function f(Ef) as 

(1-C 2 ) 2 1 



f(E f ) = ±8 ±1 



\K\Ef 



'(1-C 2 ) 2 



15 + 5(1-C 2 ) 2 — + ^(1-C 2 ) 4 
KE f 2 V 



1 



\K\Ef 
1 



{KE f f 



(12) 



(13) 



III. RESULTS AND DISCUSSION 

In FIG. [U we plot the central chemical potential of 
a non-rotating Fermi system calculated from Eq. (jTTJ) 
as a function of inverse scattering length for two differ- 
ent representative values of anharmonicity. We use N = 
2.0 x 10 6 number of atoms in the trap with 8 = 0.045. Our 
calculation shows a kink at unitarity limit due to the dis- 
continuity in the function e(y) proposed by Manini and 
Salasnich |l3[ (This kink appears in all the calculated 
macroscopic quantities). 

Solving Eq. {SI for the case of harmonic potential trap 
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FIG. 3: The lowest order axial breathing mode frequencies 
of a rotating Fermi gas in a harmonic trap with aspect ratio 
S = 0.045. The rotation frequencies are £ = (solid line), 0.5 
(long dashed line) and 0.9 (short dashed line). 

(K = 0) in a rotating (£ ^ 0) Fermi system, the lowest 
axial and radial breathing modes frequencies are given 
by 

< = |{7(2 + 5 2 - 2C 2 ) + 2(1 + 5 2 + C 2 ) 
±{[ 7 (2 + 5 2 - 2C 2 ) + 2(1 + 5 2 + C 2 )] 2 - 

85 2 [ 7 (-3 + C 2 )-2(l + C 2 )]} 1/2 }- (14) 

For the isotropic trap (5 — 1), at non-interacting limit 
and at unitarity limit (7 = 2/3), the mode frequen- 
cies are ui m = 2 and uj m = \/2 + 2Q 2 /"i. In the deep 
BEC limit (7 = 1), the mode frequencies are u> rn = 

^(1/2) (7 ± \/9~— 8C 2 )- For the case of harmonic poten- 
tial (K — 0) in a non-rotating (£ = 0) limit, Eq. (TM]) 
reduces to 

c 2 „ = i{2(l + ( 5 2 )+7(2 + 5 2 ) 

±V[2(1 + S 2 ) + 7(2 + S 2 )} 2 - 8(2 + 37)<S 2 }. (15) 

For the case of highly anisotropic limit (5 << 1), the 
two mode frequencies at 7 = 2/3 are io\jui z = \J\1jh 
and Lo-2.jui r = ydO/3 as expected. For this case, the 
two mode frequencies at 7 = 1 are uj\/lo z = y^5/2 and 
W2/w r = 2. In the BCS-BEC crossover region, the lowest 
order breathing modes frequencies are calculated from 
Eq. (fl~5|) by using the 7 from Eq. ([9]). The results for 
several representative values of Q are given in FIG. and 
FIG.H 

We solve Eq. ([8]) for the breathing mode frequen- 
cies for various values of K in both rotating and non- 
rotating Fermi systems. We calculate the central chem- 
ical potential /io for fixed number of atoms N — 
2.0 x 10 6 in the trap. We find that as d = [1/(1 - 
( 2 ) 2 }(Kh/M 2 u3)(ti Q /hiu r ) = [1/(1 - ( 2 ) 2 ]Kjl increases, 
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FIG. 4: The two lowest radial breathing mode frequencies as 
a function of d — [1/(1 — C, 2 ) 2 ]Kfl. Black (lowest mode) and 
gray (second lowest) solid lines are the mode frequencies for 
weakly interacting limit and unitarity limit (7 = 2/3). Black 
(lowest mode) and gray (second lowest) dashed lines are the 
mode frequencies for deep BEC limit( 7 = 1). The value of 
S — LUz/uJr ~ 0.045. The two modes merge at some critical 
d and beyond this critical value, the atomic cloud is unstable 
against the breathing mode oscillations. 
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FIG. 5: Axial breathing mode frequencies as a function of K 
for £ = (black), 0.3 (long dashed), 0.6 (short dashed) and, 
0.9 (dotted) at unitarity. We use the values S = 0.045 and 
N = 2.0 x 10 6 . 

the lowest order radial breathing mode frequency in- 
creases while the next lowest order breathing mode fre- 
quency decreases. Then at a critical value of d = d c , 
these two modes merge and beyond this critical d c , the 
atomic cloud is unstable against the oscillations. FIG. H] 
shows the lowest and the next lowest order radial mode 
frequencies at 7 = 2/3 (non- interacting limit and unitar- 
ity limit) and 7 = 1 (deep BEC limit). 

As evidenced by FIG.O the lowest order axial breath- 
ing mode frequencies are almost insensitive to the radial 
anharmonicities. FIG. [S] shows the axial breathing mode 
frequencies for various values of rotational frequencies as 
a function radial anharmonicities at unitarity. 

In the BEC limit where 7 = 1, we calculate the lowest 
and next lowest radial breathing mode frequencies as a 
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FIG. 6: Radial breathing mode frequencies at BEC limit 
(7 = 1) for £ = (solid line), 0.3 (long dashed line) and 
0.5 (short dashed line). The black lines are lowest breathing 
mode and the gray lines are second lowest breathing mode in 
an anharmonic trap. The value of 8 — u z /uj r — 0.045 and the 
atom number is N = 2.0 x 10 6 . 
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FIG. 7: The two lowest radial breathing mode frequencies at 
for K = 0.005 (solid line), 0.01 (long dashed line) and 0.05 
(short dashed line). The black lines are lowest breathing mode 
and the gray lines are second lowest breathing mode in an an- 
harmonic trap. The black dotted line is the lowest breathing 
mode frequency in a harmonic trap. The gray dotted lines are 
the breathing mode frequencies at K = —0.005. The value of 
5 = LJz/La r = 0.045 and the atom number is N = 2.0 x 10 6 . 

function of radial anharmonicity K for three represen- 
tative values of rotational frequencies £. We fixed the 
number of atoms to be 2.0 x 10 6 and 5 = u>z/w r = 0.045. 
As shown in FIG. as we increase K the lowest or- 
der radial breathing mode frequency increases, while the 
next lowest order radial breathing mode frequency de- 
creases. Further increase of K merges these two modes 
and beyond this merging point the atomic cloud is un- 
stable against the oscillations. 

The lowest order and next lowest order radial breath- 
ing mode frequencies as a function of the interaction pa- 
rameter [l/(fc/a)] are shown in FIG. [Jj In the presence 
of radial anharmonicity, the lowest order mode frequency 
tends to increase in the BCS regime while the next lowest 
order breathing mode frequency tends to decrease. This 
deviation becomes large as the anharmonicity increases. 
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FIG. 8: The lowest order radial breathing mode frequencies 
at for K = 0.005, 0.025 and 0.001 (Top to bottom). The 
black lines is the lowest order breathing mode in a harmonic 
trap. The value of 5 = uj z /w t = 0.045 and the atom number 
is N = 2.0 x 10 6 . The dots are the experimental data for 
comparison [l(J. 



As we have discussed before, at larger Kfio values the 
two lowest order modes merge and the cloud is unsta- 
ble against the breathing mode oscillations beyond this 
point. The data in the FIG. [8] shows the same informa- 
tion as FIG. but we plot only the lowest order radial 
breathing mode frequencies for small quartic potentials 
together with experimental data from ref. [Iol |. 



IV. CONCLUSIONS 

We have discussed the breathing mode frequencies of 
a rotating Fermi gas trapped in a harmonic plus ra- 
dial quartic potential. We find that the radial breath- 
ing mode frequencies strongly depend on the rotation 
and anharmonicity through parameter d = [1/(1 — 
(, 2 ) 2 }{KTi/M 2 ujl)(iJL /Tiu! r ). As d increases, the lowest or- 
der radial breathing mode's frequency increases and the 
next lowest order mode decreases. Beyond some criti- 
cal d c , these two modes merge and the cloud is unstable 
against the oscillations. 

As the chemical potential is large in the intermediate 
BCS regime, even with a very weak quartic potential the 
parameter d is large. As a result, the lowest order breath- 
ing mode frequency increases in the intermediate BCS 
regime. Even though the Gaussian optical trap potential 
provides a negative anharmonic term in the trapping po- 
tential, this positive anharmonic behavior has been seen 
in recent experiments 0, [l(| ■ In the weak coupling BCS 
limit, the chemical potential is even larger so that we 
find the atomic cloud is unstable against the oscillations 
at large positive anharmonicities. For negative quartic 
potentials, the breathing mode frequencies tend to de- 
crease in the BCS-BEC crossover region. 
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